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We study the collision of a brane with a black hole. Our aim is to explore the topology changing 
process of perforation of a brane. The brane is described as a field theoretical domain wall in the 
context of an axion-like model consisting of a complex scalar effective field theory with approximate 
U(l) symmetry. We simulate numerically the dynamics of the collision and illustrate the transition 
from the configuration without a hole to the pierced one with the aid of a phase diagram. The process 
of perforation is found to depend on the collisional velocity, and, contrary to our expectation, we 
observe that above a critical value of the velocity, the black hole has no chance to perforate the 
wall. That is: high energy collisions do not assist piercing. We also show that, only when the model 
parameters are fine-tuned so that the energy scale of the string is very close to that of the domain 
wall, the collision of the wall with the black hole has a possibility to provide a mechanism to erase 
domain walls, if the hole expands. However, in such cases, domain walls will form with many holes 
edged by a string and therefore disappear eventually. Therefore this mechanism is unlikely to be a 
solution to the cosmological domain wall problem, although it may cause some minor effects on the 
evolution of a domain wall network. 



I. INTRODUCTION 

A successful theory of brane collisions would provide 
an essential tool to understand many features, and, on a 
practical level, to make specific predictions, in the con- 
text of cosmology in string/M-theory scenarios, where 
the Big Bang may be associated to the collision of a pairs 
of branes. Many aspects of the problem are discussed in 
Ref. Q. 

The theory of brane collision is being developed mainly 
in two directions: a geometrical one, based on the DBI 
action and the theory of minimal surfaces, and the other 
within a super-gravity approach, by (exactly) solving the 
time-dependent equation of motion. For both approaches 
see Ref. [l[ and the list of references given there. The 
problem is also being addressed numerically within an 
effective field theory context, where the branes are mod- 
elled as domain walls (some examples are Refs. Sli). 

In all the works mentioned above, the collision takes 
place between two branes immersed in a higher dimen- 
sional space. Here we would like to modify the starting 
point, and ask what happens when a brane collides with 
another 'object', that is not necessarily a brane. The 
object we have in mind is a black hole. 

From the point of view of the brane-uni verse, this is 
a natural question to ask, and physical set-ups of this 
sort can be easily imagined at early stages of the evo- 
lution of the universe, where perturbations might have 
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collapsed and formed primordial black holes. Besides, 
the study of this problem can be motivated from var- 
ious other directions. A recent one lies in the context 
of topology changing transition in gravitational systems. 
A well studied example of this sort is provided by the 
{merger) transition between a Kaluza-Klein black hole 
and a non-uniform black string [5, 6]. Such an example 
presents striking similarities with the system of a brane 
and a black hole, as it has been recently highlighted in 
Ref. Q, where the importance of exploring these phe- 
nomena in the context of field theory has been stressed. 
The reason is that when a change in the topology of the 
spacetime occurs, the curvature is expected to diverge. 
Thus, a field theory model would provide a reasonable 
example to explore the process of topology change in a 
physical situation free of singular behaviors. Reference 
[7| proposes a working-model to explore these similari- 
ties: a (Dirac-Nambu-Goto) brane attached to a black 
hole: a brane-black hole system. Related study has also 
been carried out by the present authors both within a 
geometrical approach [8|1 , and in the context of domain 
walls in field theory [jjl llOl] . 

In this work we wish to explore some other features of 
the brane-black hole systems, which are important when 
field theory constructions are employed: when the brane 
is described as a field theoretical domain wall, topology 
changing processes depend on the specific field theory 
model used. Here we will show this explicitly and con- 
sider an example where a 'topological process', which is 
impossible in the set-up considered in Refs. @, d, B EH , 
can be realized: the perforation of a brane due to the 
collision with a black hole and subsequent nucleation of 
a string loop on the brane world-sheet. This is evidently 
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relevant in the context of the brane universe, in the terms 
described at the beginning of this section. For this reason 
we consider the problem in general dimensions and dis- 
play the results of numerical simulations in 4D, 5D and 
6D. 

Let us now turn to the problem at hand: the perfora- 
tion of a brane. A similar question was investigated in 
flat space in Refs. [HI, [ll]- The main point was to show 
that when there is an approximate gauge symmetry or 
a series of phase transitions, hybrid topological struc- 
tures may arise. One example of this sort is provided 
by domain walls attached to strings, and concrete real- 
izations are provided by axion models, having an exact, 
spontaneously broken discrete symmetry, embedded 
in a t/(l) group. The vacuum structure presents, due 
to the discrete symmetry, a degeneracy, and an axion 
domain wall is a minimum energy configuration which 
interpolates between two neighboring vacua. In some 
cases (N = 1), the walls are topologically unstable by 
quantum nucleation of a hole, whereas N > 1 models 
are, instead, stable. Despite of this difference, the decay 
rate of N = 1 walls is exponentially suppressed. Further 
work, carried out in Ref. [l3j], considered gravitational 
effects and has shown that, as far as 'de Sitter' walls are 
concerned, the nucleation of one hole is not sufficient to 
completely destroy the wall, and at least four holes are 
necessary for this purpose. Interestingly, as the authors 
of Ref. [13j suggest, this may signal a more generic kind 
of instabilities of p-branes in supergravity models. 

In the present article, we wish to discuss the formation 
of a puncture on the domain wall as a classical phenom- 
ena, occurring due to the collision with a black hole and 
discuss the topological stability from a classical stand- 
point. After recalling the main features of the field the- 
ory model we employ, we will report on our results about 
the dynamics of the collision. The complex scalar effec- 
tive field theory of Ref. Qjl is used to describe the field 
theoretical domain wall. Such model, in principle, allows 
one to consider the phenomena of perforation. The col- 
lision between an axion domain wall and a black hole is 
then simulated numerically illustrated with the aid of a 
phase diagram (Fig. [3]) . 

Before jumping into the calculations we would like to 
comment on another cosmological scenario for which the 
process we are about to study is important: the cosmo- 
logical domain wall problem. It is not our aim to discuss 
it in detail (a thorough discussion and a list of references 
can be found in Ref. [3), rather, we would like to make 
some comments on a conjectured solution, proposed in 
Ref. [H. 

It is well known that domain walls occur in any gauge 
theory in which a discrete symmetry, not part of the 
gauge symmetry, is spontaneously broken. Hence, it is 
difficult to imagine that domain walls, regardless of what 
the fundamental theory may be, have not been created in 
the early universe, their presence has, in most cosmolog- 
ical models, disastrous consequences. This is, essentially, 
due to the fact that domain walls carry too much energy 



and destroy the prediction of standard Big Bang model 
(T^ | . The question is how to avoid their domination over 
the energy density of the universe. 

This problem has been explored extensively, and a 
range of possible solutions has been considered. Inflation 
is a popular way out, since it can dilute, over an expand- 
ing volume, the energy density of the walls. However, this 
mechanism may not work if the walls are created at en- 
ergies below the inflationary scale typically of order 10 16 
GeV. Other remedies have been suggested. One example 
is the 'bias mechanism', discussed in Refs. [19|, 12CJ, l2l| . 
which introduces a small difference in the energy density 
of the vacua on the two sides of the wall. As a result, 
the wall may collapse. This mechanism, however, suffers 
from some limitations, and in realistic models, instead of 
decaying, the walls may form stable bound state struc- 
tures (22J. High temperature symmetry non restoration, 
dissipation due to interaction of the walls with the sur- 
rounding matter, and many other possibilities have been 
discussed in the literature (see [18[ for an extensive re- 
view), but, to the best of our knowledge, none of the 
possible remedies proposed so far seem to have a 'uni- 
versal' character that cures the domain wall problem in 
a general and model independent way. This leaves space 
for discussion of other mechanisms, and here we would 
like to comment on the one proposed in Ref. 17]. The 
idea is that, interacting with primordial black holes, the 
domain walls are perforated and a hole forms. Since the 
process is classical, and thus may not be suppressed, the 
subsequent expansion of the hole may provide, in princi- 
ple, an efficient way to get rid of the walls. 

Although, at first glance, the mechanism may seem at- 
tractive, when considering specific realizations, one can 
easily argue that it may only work in cosmologically rel- 
atively safe situations. The reason is simply that a hole, 
when created, must be bounded by a string. In other 
words, we must consider an hybrid domain wall-string 
system, like those arising in axion models, mentioned 
above. Evidently, the mechanism does not work in the 
case of N > 1 (N essentially tells us how many 'do- 
main wall branches' are attached to the string), since the 
production of a hole is topologically impossible, i.e. in 
the N > 1 case the different domains are in different 
vacua. If a hole is formed, one can go from one vacuum 
to another without going over a potential barrier, which 
is impossible for N > 1. In the case of A = I, there 
is only one vacua, and the interpolating field configura- 
tion goes around the bottom of the potential. The meta- 
stability of such walls makes them harmless for cosmol- 
ogy, since the walls rapidly break up as a result of multi- 
ple self-intersections and decay into elementary particle 
and gravitational waves, as discussed in Ref. [l2|. A re- 
cent review on axion cosmology can be found on Ref. [23| . 

The previous considerations make the mechanism pro- 
posed in Ref. [13] less appealing (one may also try to 
imagine more complicated topological processes in the 
case N > 1 , however the simulation we carried out in the 
case N — 2,3 suggest the only recombination will take 
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place). However, if the energy scale at which strings form 
is much larger than that at which domain walls form, in- 
flation can occur in between, and usual dynamical pro- 
cesses will not be an efficient way of destroying the walls. 
It is, thus, reasonable to ask whether the collision with 
primordial black holes may enhance the destruction of 
the wall and compete with other mechanisms in specific 
examples. 

Anticipating the results we will discuss in this paper, 
the mechanism turns out to be efficient only if the do- 
main walls form at energy scales immediately below those 
at which strings form. However, in this case the domain 
walls have a lot of holes edged by a sting from the begin- 
ning and therefore are cosmologically harmless. For this 
reason the mechanism of perforation turns out to be less 
attractive than it might seem at first sight. 



II. AXIONIC BRANES 

Axionic domain walls were introduced after Peccei and 
Quinn proposed a solution to the strong CP problem [lj| , 
and were studied in the context of topological defects and 
cosmology (See for example Refs. [HI, E3, EE The 
Peccei-Quinn symmetry is a global U(l) symmetry which 
is broken at an energy scale fg typically constrained by 
accelerator experiments and astrophysics, leaving the en- 
ergy scale fg to vary within the so called axion window, 
between 10 10 GeV and 10 12 GeV. This corresponds to 
constraining the axion mass roughly within the range be- 
tween 10 -5 eV and 10 -2 eV. Due to the U(l) symmetry, 
global strings are produced during the symmetry break- 
ing. The axion, initially massless, becomes massive at 
the QCD scale due to instanton effects. When the tem- 
perature, after reheating, reaches values of order of 100 
MeV, the potential acquires a non vanishing term: 



U N = CmJ(l- cos N0/fe) 



(1) 



with 9 being the axion field, C a constant of O(l) and N 
an integer which depends on the detailed structure of the 
model. The above potential has a Zn symmetry which 
accounts for domain wall solutions. 

Here we will consider a (p + 2)— dimensional effective 
theory described by the following action: 



S = / a? +2 x L 



where 



where 



(2) 



(3) 



U(§, *t) = j ($$t -rj 2 ) 2 - 2/r 2 $ 2 (cos 9-1). (4) 

The above effective theory, corresponding to the axion 
model with TV = 1, will be our working model. The 



complex scalar field $ is parametrised as $ = p e and 
9 represents the axion field. This model has a unique 
vacuum state at |$| = ry, 9 — and no exact symmetry. 
However the model has an approximate U(l) symmetry 
for ii 2 -C A?7 2 , and there are accompanying global string 
solutions. If the symmetry were exact, the phase 9 would 
change uniformly around the string, but the second term 
in the potential forces 9 to zero, so that all the variation 
of 9 from to 2-7T is confined to a domain wall centered 
at 9 = it, whose boundary is attached to a string. The 
thickness of the domain wall ~ is much larger than 
that of the string core ~ (\f\rj)^ 1 . To see the presence of 
domain wall solutions, it is sufficient to notice that away 
from the string core, the VEV of $ is < $ >= fye 10 , and 
the above Lagrangian becomes 

L = {-d^d^O - 2^ 2 (cos(9 - I)) rj 2 , 

which is of the form of a sine-Gordon model. In such 
a case, plane-symmetric domain wall solutions are well 
known to be of the form 

9{x) = 4arctane AIZ , 

where z is a coordinate perpendicular to the wall. 

We record here some useful scaling properties of the 
above N = 1 axion model. Notice that upon a coordinate 
dilatation, 



followed by 



and 



the action scales as 



o/a , 



a rj \ 

2 2 

a ri 



S ^ rr~ S . 

To 



(5) 



It is trivial to see that: we are free to rescale r] to unity 
(in four dimensions, energies will be measured, then, in 
fractions of rf); and that a coordinate dilatation at the 
classical level corresponds to simultaneous variation of 
the coefficients A and /i 2 . 

We conclude the section by commenting on the gen- 
erality of the model considered here. It can be easily 
understood that for the process of perforation to occur, 
the essential feature is that the brane must be metastable, 
meaning that the formation of a hole bounded by a string 
has to be allowed. Therefore the potential must be char- 
acterized by at least three quantities: one giving the en- 
ergy scale at which the domain wall forms; the second 
being the size of the vacuum manifold that accounts for 
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string formation; and, finally, the height of the potential 
leading to string formation. All these features are com- 
prised in the above model j4}, which depends on three 
(free) parameters (77, A and fi) that completely charac- 
terize the above three quantities. 

The above scaling properties turns out to be a conve- 
nient bonus when considering the interaction between a 
black hole and the domain wall. Defining the ratio of the 
portion of the wall accreted by the black hole and the 
mass of the black hole as 

<T W R P 2 r, 

e — = ii-n R , 

Rp -i P / ' 

the probe brane approximation can be written as 

e«l . 

Here R is a typical scale of order of the black hole horizon 
and o w ~ /wy 2 is the tension of the wall. Now, in view of 
the above scaling properties, the quantity e changes as 

e -> e l r ll ■ 

It is then immediate to understand that the above scaling 
essentially enforces the probe brane approximation. 

It is also easy to see that in a cosmological context 
the probe brane approximation is valid. If we assume 
that the wall do not dominate the energy density of the 
universe, then the following inequality follows (from the 
Friedmann equation): 



with 



< 1 



Assuming that the size black hole is much smaller the 
Hubble radius, one recovers the relation e <C 1. Having 
said this, we can now safely move on to the next section. 



f=l- 



r p- 



(6) 



Here we could set the horizon radius th to unity without 
loss of generality because the scaling property mentioned 
above. Then the Klein-Gordon equation takes the form 

0= - /-V± + ^ r (r'/W±) 

+ , > (shrP- 1 ^±)+^± , (7) 



where we have defined $ = "0+ + iip- and 



U+ - iU- 



dU_ 



(8) 



The collision process, in the probe brane approximation, 
is then described by the above time-dependent system of 
partial differential equations with appropriate boundary 
and initial conditions. 

Initially the domain wall is located at a distance zq 
from the equatorial plane of the black hole. The initial 
configuration can be obtained by using, for example, the 
relaxation method. This solution will not differ much, as 
intuition suggests, from the static solution in flat space 
if zq is large enough, i.e., if the brane is far away from 
the black hole. And, in fact, the relaxed solution, can be 
compared with the 'flat space' configuration in the limit 
/i 2 < Xt] 2 : 



p(t = 0,r,6) = V , 6(t = 0,r, 



4 arctan e 



Here z = rcos8. The motion of the domain wall at the 
initial stages of its evolution (still far away from the black 
hole) is described by boosting the solution 



III. COLLISIONS WITH BLACK HOLES 



^i ) ( 1 ((z-z )-vt)) 



Axionic branes are metastable, in the sense that a 
string loop can be nucleated spontaneously and, in princi- 
ple, can destroy, at least locally, the wall. However, since 
the nucleation rate is exponentially suppressed (and even 
if a hole is created, the brane may inflate with a higher 
expansion rate than the hole), it is not clear whether 
the brane survives on cosmological time-scales. For this 
reason, we want to investigate a classical (and thus not 
suppressed) mechanism of hole nucleation, as the colli- 
sion with a primordial black hole may be. This process 
is described by the Klein-Gordon equation for the com- 
plex scalar field $, 



dU 







We consider the spacetime of a (p + 2)-dimensional black 
hole, 



ds z = -fdi z + f^dr 2 + r 2 dn 2 p , 



where 7=1/ Vl — v 2 is the Lorentz factor. The previous 
expression will describe the initial conditions. In order to 
determine the evolution of the system we also need to fix 
the boundary conditions, that will describe the motion 
of the black hole during the collision far away from the 
hole. For the numerical simulation 'far away' means sev- 
eral times the horizon size. In the numerical simulation 
we will take a relaxed flat space solution ip^\z, t) as a 
configuration to describe the initial and infinity bound- 
ary conditions. At the horizon the boundary conditions 
are written as 

dtil>±(t,r H ,6) = , 

since the lapse vanishes there in our coordinates. 

To study the evolution of the system we have to solve 
the system of PDE ([7]) . To impose the horizon boundary 
conditions, spherical coordinates are appropriate. Then 
we need to treat near the z-axis a little carefully. A conve- 
nient way to handle this problem is to use the harmonic 
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decomposition for the angular coordinate. We initially 
decompose the solution in terms of a complete basis of 
smooth global functions in the 9 direction: 

n 

i>±(t,r,6) = £u±(i,r)^ fe (cos0), 

k=0 

where n is the number of harmonics used. We choose the 
basis functions (f>k(cos9) to be the Chebyshev polynomi- 
als of the first kind, which are regular at 9 = 0, n and 
satisfy the following orthogonality relation (x — cos 9) : 




<j>i(x)<j>j{x)w{x)dx = —Sij , 



with w(x) — l/v'l — x 2 . (Notice that the normalization 
for 4>o differs from the conventional one.) Using the previ- 
ous decomposition, it is possible to rewrite the evolution 
equations (|7|) as follows: 

4 = (r*fd r uf) + ^J Jh ■ «± , 0) 

where 

Iij = J w(x) [d e (f>i(x) + (p - 1) cot 9de<pi (x)] <j)j(x)dx , 
and 

r+l 

Vf = J U±<j> 5 {x)w{x)dx . 

The last step of the numerical procedure that we adopt 
consists of solving the previous equation @ by means 
of standard finite difference methods. As a check on our 
numerical scheme we have evolved the domain wall in 
flat space, modifying appropriately the boundary condi- 
tion at the horizon, making sure that the wall propagates 
smoothly, as it must. 

IV. DISCUSSION OF THE NUMERICAL 
RESULTS 

Before turning to the numerical results, it is useful to 
discuss what is the physically relevant range of parame- 
ters. First of all, it can be easily argued that for values 
of A below a critical value, domain wall solutions cease 
to exist. To see this, it is sufficient to notice that the 
potential has an approximately U(l) symmetric string of 
vacua. Along the trajectory defined by 5R($) = ip and 
3(<I>) = (see Fig. [T]), the potential assumes the form 

where 9{x) is the step function. The previous function 
has two minima, one at <j>T — —7], the 'true' one, and one 

at 4>f = a/?/ 2 — 8^-, the 'false' one. Roughly speaking 



the domain wall configuration interpolates between these 
two vacua. It is then easy to show that the minima at <j>p 
becomes a saddle point for A ~ 8[i 2 /t] 2 . In this case the 
approximate U(l) symmetry is absent and it is obvious 
that domain wall solutions do not exist. Once the values 
of rj and [A are fixed, 8fj?/r] 2 gives an upper bound on 
the critical value Ai for A below which there is no static 
stable domain wall solution in flat spacetime. 



D 




FIG. 1: Potential in the N = 1 case. 

Here it is worth discussing an alternative way to esti- 
mate the critical value Ai. Let us first focus on the flat 
space case. When we consider a configuration with a hole 
bounded by a string, the hole may expand or collapse, 
depending on its size. If there is a stable domain wall 
configuration without a hole, we can argue that, if a hole 
smaller than a critical size forms, then it will collapse. 
On the other hand, if the hole forms with a size larger 
than the critical one, it will expand. Hence, between 
these two possibilities, there must be a marginal (unsta- 
ble) configuration at which the hole neither expands nor 
collapses. Conversely, if there is no unstable static con- 
figuration with a hole bounded by a string, the domain 
wall is connected to the expanding hole without a bar- 
rier. In this case the domain wall configuration without 
a hole will not be stable. 

To see if there exists an unstable static configuration, 
we make use of the thin wall (Dirac-Nambu-Goto) ap- 
proximation. In this approximation, we will find that 
this critical configuration always exists. However, once 
we take into account the finite thickness of the string, 
such a solution will be expected to vanish when the sting 
thickness becomes as large as the radius of the hole. The 
thickness of the string r s can be estimated from the cur- 
vature of the scalar field potential as r s rs a/\/Xr}, where 
a is a constant of 0(1). To evaluate the radius of the hole 
for the critical configuration, we extremize the energy of 
the the static configurations consisting of a pierced brane 
and a string bounding the hole with radius r, 

/>oo 

E cx n s r p - 1 +a w drr^ 1 , (10) 

J r 

where fi s as rj 2 and a w w firf 1 are the string and the wall 
tensions, respectively. The variation of the energy gives 
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the radius of the hole at the critical configuration as 



Thole 



r = (p- 1) 



Ms 



t 

1 



(11) 



where we have introduced a constant f3 of O(l). Let us 
assume that the condition for a critical configuration to 
exist, after taking into account the string thickness, is 
given by r s < KThoie with a constant k of 0(1). Then, 
the threshold value for the stability of an un-pierced wall 
is estimated as 



Ai 



9 2 

or pr 
P 2 k 2 rf 



(12) 



The dependence on fi and r\ is the same as before, ro 
marks a critical values, above which a hole will expand 
and below which it collapses. 

The above analysis has the advantage that it can be 
easily extended to the black hole background, which we 
will now consider taking the domain wall on the equato- 
rial plane of the black hole. 

Even if the scales at which strings and domain walls 
form are very different, and usual dynamical mechanisms 
do not help to destroy the walls, collisions with black 
holes may still have a chance to be an effective mecha- 
nism, because, when the horizon radius of the black hole 
is larger than the size of the critical hole, ro ss j3/n, the 
domain wall will be perforated. Cosmologically we can 
expect this to be reasonable if formation of primordial 
black holes is delayed and occurs at a low scale. This 
suggests that whether the domain wall is pierced or not 
seems to be more or less independent of A. However, this 
naive expectation is wrong. To show this, we repeat the 
analysis that we have done above under the presence of 
a black hole. The expression for the energy (fTO]) is to be 
replaced with 



By varying this expression with respect to r, we find 



0=1- -r 1 ^ 
2 



1 - 



1 



r p i rQ 



The coordinate radius of the hole Thole is determined by 
solving this equation. 

However, rather than Thole, a more appropriate quan- 
tity to be compared with the string thickness r s is the 
proper distance between the rim of the hole and the hori- 



A(r ) = 



f- 1/2 (T)dr 



We can reasonably imagine that, once the string thick- 
ness r s becomes larger than A, the static configuration 
of a brane with a floating string will not exist. Then, 
following the same discussion as in the case of the flat 



background, domain wall configurations with a hole of 
any size will become unstable. Comparing r s with A(ro), 
it is easy to see that this occurs for 



A > A 2 



a (i 



(3 2 k 2 t] 2 A 2 (r )' 



This tells us that there exists a threshold value A2, above 
which the black hole cannot pierce the wall. We find 
that A2 differs from Ai by the factor r\ / A 2 (ro) , which is 
plotted in Fig. [5] against logr . Independent of the value 
of ro, this ratio does not deviate much from unity. When 
ro > 1 (recall that the horizon radius is normalized to 
unity), the existence of a black hole can be neglected. 
Then the ratio becomes unity and we have Ai « A 2 . 
On the other hand, in the limit ro <C 1, we have A2 ~ 
(p-l)Ai. 



3 

2.5 
2 

1.5 

J 



-4 -3 



-Log r D 



FIG. 2: rg/A 2 vs logr . For small ro, increasing values of the 
dimension, increases the ratio r^/A 2 . The curves refer (from 
bottom to top) to D — 4, 5, 6 spacetime dimensions. 

In summary, according to the above arguments, we can 
expect that for A < Ai, there will be no static domain 
wall configuration; for A > A2 the domain wall configu- 
ration exists and cannot be perforated in the adiabatic 
process. In the intermediate range Ai < A < A 2 the do- 
main wall exists at a distance from the black hole and 
can be perforated eventually as it gets closer to the black 
hole. 

It is now time to present the numerical results. How 
our simulation goes is summarised in Figs.OHl where var- 
ious snapshots of the evolution of the domain wall dur- 
ing the collision are illustrated. We present both cases 
where a hole is created and expands or only reconnec- 
tion takes place. The figures show the evolution of the 
field components and of the energy density, and, in some 
case, to emphasize the formation of the hole, we have 
presented the plots of the energy density from two differ- 
ent viewpoints. The parameters of the simulation have 
been clipped to n — 40 harmonics, 400 grid points and a 
maximal radial distance up to 20 times the horizon size. 
To scan the parameter space we fixed the value of the 
symmetry breaking scale r\ to unity and varied /i and A. 
The collisional velocity has been sampled from 0.9 down 
to 0.001 in units of the velocity of light. 

Focusing on whether a hole is formed or not, the re- 
sults of simulations for various choices of parameters are 
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summarized in Fig. [3l The dots in Fig. [3] are calculated 
performing the simulation and increasing the collisional 
velocity at steps of 0.05. The critical value is chosen as 
the mean point between the last value for which a hole 
forms and the first for which it does not. The error made 
in this estimate is taken care of by drawing the dots with 
a radius of 0.025. The plots are made for three differ- 
ent values of /i; p = 0.2,0.4 and 0.8. Interestingly, and 
contrary to our intuition, in a certain region of A, we 
record the presence of a critical velocity, above which the 
formation of a hole is not observed. Quadratic best fit 
curve for the critical velocity is also shown for each value 
of /i. From this plot we can distinguish the parameter 
space into three regions, as schematically described in 
Fig. 21 (1) domain wall solutions do not exist even in 
the flat background, (2) the domain wall is pierced by 
a black hole after collision, and (3) the domain wall is 
not pierced. Above the 'critical line' the domain wall is 
stable in the sense that the black hole is unable to pierce 
the domain wall. 

The low velocity limit in this plot will correspond to 
the adiabatic case that we have discussed above. Hence, 
A2 can be read from this limit of the curves for the crit- 
ical velocity. When we prepare the initial configurations 
for the present simulations, the threshold values Ai, at 
which a domain wall solution marginally exists in the flat 
background, are obtained simultaneously. Thus the ob- 
tained numerical values of A1//1 2 are 36.25, 43.75, 45.31 
for jU = 0.2, 0.4 and 0.8, respectively. These values are al- 
most the same as the values read from the large velocity 
limit of the curves for the critical velocity. For [i = 0.2, 
ro is sufficiently large and the ratio A2/A1 is close to 
unity. While, for /z = 0.8, ro is close to unity and the 
ratio A2/A1 is close to 2. These results perfectly agree 
with the picture presented above. We have performed 
the same check in four and six dimensions, finding good 
agreement. 





domain walls cannot 
be perforated 



FIG. 4: Schematic phase diagram representing the depen- 
dence of the critical velocity from A. 



We briefly discuss the potential relevance for the cos- 
mological domain wall problem. First of all, we notice 
that, contrary to our intuition, the collision with a black 
hole can be an efficient way of getting rid of the domain 
walls in the region of parameters where the energy scale 
for the string is very close to that for the domain wall. 
This is an unexpected constraint found in the present pa- 
per. One may think that in a concrete model, the mech- 
anism may compete with other processes in erasing the 
wall. However, it is easy to see that this is not the case. 
In fact, it can be argued that generically the perforation 
mechanism is not a likely solution anyway. Cosmologi- 
cally N = 1 axion-domain walls are harmless if they are 
associated with many strings. When two energy scales for 
the string and the domain wall are very close, it is almost 
impossible to imagine that only domain walls are formed 
without forming strings. In simple words, the barrier for 
the string formation is not very large compared to the 
energy 'stored' in the walls, and it is reasonable to be- 
lieve that any dynamical process like wall-wall collisions 
or fluctuations of the wall will lead to the creation of a 
hole in the wall. This suggests that the collision with the 
black holes will not play a major role in the fate of such 
domain walls, although it may produce minor modifica- 
tions in the evolution of a domain wall - string system. 



V. CONCLUSIONS 



FIG. 3: This figure illustrate the dependence of the criti- 
cal velocity, above which a hole does not form, on the ratio 
Xrj 2 1 fj 2 . It refers to D = 5 dimensions and the curves, from 
left to right correspond to the small, intermediate and large 
fi range (fi — 0.2, 0.4 , 0.8 respectively). The points are com- 
puted from the numerical analysis by fixing the values of [i, 
and performing the simulation varying A and the collisional 
velocity. 



Since it was shown that gravity can be localized on 
a (p + 1)— dimensional submanifold of a higher dimen- 
sional spacetime, the idea that 'our' universe could be 
a 3— brane immersed in a higher dimensional bulk space 
has been a subject of intensive study. 

In this context, it is interesting, not only from an aca- 
demic perspective, to ask what happens when a collision 
with a brane occurs. One may imagine that the collision 
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occurs between two branes, as in some stringy cosmolog- 
ical models, or between a brane and a black hole, as it 
could happen, when considering primordial black holes. 
This last problem happens to be relevant also in four 
dimensional standard cosmology, when the brane is in- 
terpreted as a cosmological domain wall. 

The aim of this paper was to study such a prob- 
lem, namely the collision between a brane and a black 
hole in a topologically non-trivial model. The meaning 
of 'topologically non-trivial' is that processes, different 
from recombination, can, in principle, occur. Specifi- 
cally we concentrate our attention to the possibility that 
the brane is perforated due to the collision. This can oc- 
cur when the eventual hole is bounded by a string, and, 
as working model, we consider axion type effective field 
theories, which account for hybrid defects (domain walls 
attached to strings), that, in our treatment, will model 
the brane. 

The problem of the collision is investigated by numer- 
ically solving the Klein-Gordon equations for the brane 
in the black hole spacetime. We have scanned the pro- 
cess by varying the parameters of the potential, \x and 
A, and this also effectively takes account of varying the 
size of the black hole, due to the scaling properties of the 
Lagrangian. 

The results of the simulation are characterized, for 
fixed fi, by two critical values of the parameter A: Ai 
and A2. The first one gives a lower bound on the exis- 
tence of a domain wall solution, whereas, the other value, 
A2, marks the threshold above which the domain wall 
is stable in the sense that cannot be perforated by the 
black hole. The critical value A2 separates the case where 
an unstable static domain wall configuration with a hole 
edged by a sting on the equatorial plane of the black 
hole exists and the case in which such a configuration 
does not. When the string thickness r s exceeds A, the 
proper distance between the horizon and the string core, 
the unstable configuration does not exist. The hole can- 
not shrink and the domain wall is necessarily pierced. 
On the other hand, for r s < A a static configuration of a 
domain wall with a floating string exists. Consequently, 
the hole when created, is accreted by the black hole. As 



a result the wall is not pierced. 

In the intermediate region between Ai and A2, where 
the domain wall can be perforated, we observe the ex- 
istence of a critical collisional velocity, v cr u, (unexpect- 
edly) above which the domain wall is also stable. This 
fact is illustrated with the aid of the phase diagram of 
Fig. [3] We have numerically evaluated the ratio A2/ Ai , 
and compared it with the one estimated from analytic 
arguments using the static configuration, finding agree- 
ment. 

Finally, we have discussed the relevance of the above 
process in the context of a recently conjectured solution 
to the cosmological domain wall problem. We argue that 
the mechanism is efficient only when domain walls form 
immediately after strings. However, in such cases, string 
formation cannot be suppressed when domain walls are 
formed, because the potential barrier to overcome in or- 
der to have string production is comparable to the energy 
stored in the domain walls. Hence, those domain walls 
are quickly eaten up by the expanding holes (formed be- 
fore the collisions) and disappear anyway. Therefore it 
seems unlikely that this mechanism plays an important 
role as a solution to the the cosmological domain wall 
problem, even though may have a small effect in the evo- 
lution of domain wall networks. 
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FIG. 5: Snapshots of the evolution of ifr±. The values used in this simulation are A = 10, fj, = 0.4 and rj — l.The velocity if 
chosen to be 0.1. The black hole horizon is a circle of radius 1 centered at the 0. The spacetime dimensionality is 5. One can 
observe the formation of a hole that locally destroys the wall. 




FIG. 6: Snapshots of the evolution of the energy density p. The values used in this simulation are the ones of the previous 
figure. To emphasize the expansion of the hole, we illustrate the evolution of p from two viewpoints. 




FIG. 7: Snapshots of the evolution of tp± and of the energy density p. The values used in this simulation are A = 10, /i = 0.4 
and r\ — 1. The velocity if chosen to be 0.5. The black hole horizon is a circle of radius 1 centered at the 0. The spacetime 
dimensionality is 5. In this case only recombination occurs and the domain wall reconnects behind the black hole. 




FIG. 8: Snapshots of the evolution of tp±. The values used in this simulation are A = 10, fj, = 0.4 and rj — l.The velocity if 
chosen to be 0.1. The black hole horizon is a circle of radius 1 centered at the 0. The spacetime dimensionality is 4. One can 
observe the formation of a hole that locally destroys the wall. 




FIG. 9: Snapshots of the evolution of the energy density p. The values used in this simulation are the ones of the previous 
figure. 



